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Abstract 



The gonihedric model of random surfaces on a 3d Euclidean lattice has equivalent rep- 
resentation in terms of transfer matrix K(Q i} Qf) which describes the propagation of loops 
Q. We extend the previous construction of loop transfer matrix to the case of nonzero 
self-intersection coupling constant k. We introduce loop generalization of Fourier transfor- 
mation which allows to diagonalize transfer matrices depending on symmetric difference of 
loops and express all eigenvalues of 3c? loop transfer matrix through the correlation func- 
tions of the corresponding 2d statistical system. The loop Fourier transformation allows 
to carry out analogy with quantum mechanics of point particles, to introduce conjugate 
loop momentum P and to define loop quantum mechanics. 

We also consider transfer matrix on 4c? lattice which describes propagation of meme- 
branes. This transfer matrix can also be diagonalized by using generalized Fourier trans- 
formation, and all its eigenvalues are equal to the correlation functions of the correspond- 
ing 3c? statistical system. Particularly the free energy of the 4c? membrane system is equal 
to the free energy of 3c? gonihedric system of loops and is equal to the free energy of 2d 
Ising model. 



1 Introduction 



Various models of random surfaces built out of triangles embedded into continuous space 
R d and surfaces built out of plaquettes embedded into Euclidean lattice Z d have been 
considered in the literature [1], 0, |15| . These models are based on area action and suffer 
the problem of non-scaling behaviour of the string tension and the dominance of branched 
polymers ||. Several studies have analyzed the physical effects produced by rigidity of the 
surface introduced by adding dimension-less extrinsic curvature term to the area action 

1 

1 

T ■ S(area) H — F( extrinsic curvature), 
a 

where T is a string tension, a is dimension-less coupling constant and F(ext curv) oc 1 is 
dimension-less functional. Comprehensive review of work in this area up to 1997 can be 
found in 0. 

In || § authors suggested so-called gonihedric model of random surfaces which is also 
based on the concept of extrinsic curvature, but it differs in two essential points from the 
models considered in the previous studies. First it claims that only extrinsic curvature 
term should be considered as the fundamental action of the theory 

S = m ■ A(extrinsic curvature), 

here m has dimension of mass and there is no area term in the action. Secondly it is 
required that the dependence on the extrinsic curvature should be such that the action 
will have dimension of length A(extrinsic curvature) oc length. This means that 
the action should measure the surfaces in terms of length and it should be proportional 
to the linear size of the surface, as it was for the path integral. The last property will 
guarantee that in the limit when the surface degenerates into a single world line the 
functional integral over surfaces will naturally transform into the Feynman path integral 
for point-like relativistic particle (see Figure 1,2) [] 



m A xy (extrinsic curvature) 




As it was demonstrated in || quantum fluctuations generate the area term in the effective 
action 

S e ff = m A{extrinsic curvature) + T e ff S(area) + 

that is nonzero string tension T e ff. Therefore at the tree level the theory describes 
free quarks with string tension equal to zero, but now quantum fluctuations generate 
nonzero string tension and, as a result, quark confinement ||. The theory may consistently 
describe asymptotic freedom and confinement as it is expected to be the case in QCD. 

Our aim is to find out what the continuum limit of this regularized theory of random 
surfaces is, and if indeed the continuum limit corresponds to nontrivial string theory, does 
it describe QCD phenomena. This is a complicated problem and one of the possible ways 
to handle the problem is to go further and to formulate the theory also on a lattice. An 
equivalent representation of the model on Euclidean lattice Z d has been formulated in |7| . 
This representation depends on dimension d of the embedding space. In three dimensions 

lr This is in contrast with the previous proposals when the extrinsic curvature term is a dimension-less 
functional F (extrinsic curvature) oc 1 and can not provide this property. 
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Figure 1: It is required that the action A xy should measure the surfaces in terms of 
length, as it was for the path integral. When a surface degenerates into a single world 
line we have natural transition to a path integral . 

it is a spin model in which the interaction between spins is organized in a very specific 
way (see equation (|8]) ).In high dimensions it is a gauge spin system which contains one, 
two and three plaquette interaction terms. 

In this article I shall consider the above model of random surfaces embedded into 
3d Euclidean lattice Z 3 . The reason to focus on that particular case is motivated by 
the fact that one can geometrically construct the corresponding transfer matrix |12[ and 



find its exact spectrum (T3|. The spectrum of the transfer matrix which depends only 
on symmetric difference of initial and final loops has been evaluated exactly in terms of 
correlation functions of the 2d Ising model in [T3| . This transfer matrix has the form W2 



K K=0 {Qi, Q 2 ) = exp{-(3 [k(Q!) + 2J(Qi A Q 2 ) + k{Q 2 )} }, 

where Qi and Q 2 are closed polygon-loops on a two-dimensional lattice, k(Q) is the 
curvature and l(Q) is the length of the polygon-loop Q Q. This transfer matrix describes 
the propagation of the initial loop Qi to the final loop Q 2 . 

Our aim in this article is to extend these results in several directions. First to construct 
the transfer matrix for more general case of nonzero self-intersection coupling constant k 

K{Q t , Q 2 ) = K K=0 {Q 1 , Q 2 ) ■ exp{-4K(3 [k int (Q x ) + 2J(Qi n Q 2 ) + k int {Q 2 )} }, 

where ki nt (Q) is number of self- intersection vertices of the loop Q. For k = this expres- 
sion coincides with the previous one. In the limit k —>■ oo, one can see that on every time 
slice there will be no loops with self-intersections and we have propagation of nonsingular 
self-avoiding oriented loops. 

2 We shall use the word "loop" for the " polygon- loop" . 
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Secondly, I will demonstrate that diagonalization of any transfer matrices K(Qx AQ2), 
which depend only on symmetric difference of loops Qi A Q 2 can be performed by using 
generalization of Fourier transformation in loop space. This transformation will be defined 
by using the loop eigenfunctions 

which are numbered by the loop momentum P and s(Q) is the area of the region with 
boundary loop Q. They are in a pure analogy with plane waves in quantum mechanics 
i/j p (x) = e vpx . We shall prove that any transfer matrix K(Qi A Q2) can be diagonalized 
by using this loop Fourier transformation 

K(Pi, Pz) = E K (Qi A Q2) e^ p ' nQl) - ins ^ nQ2) = A Pl 6{P U P 2 ), 

here Ap are eigenvalues depending on loop-momentum P. Because the transfer matrix 
plays the role of the evolution operator, one can define a corresponding Hamiltonian 
operator as H(P) = —InA(P). At this point one can talk about loop quantum mechanics 
which should be defined by using conjugate loop operators Q and P and the Hamiltonian 
H(P-.Q). 

I shall also consider the model of random three-dimensional manifolds which are build 
by gluing together tetrahedra in continuous space R A and out of 3d cubes in 4d Euclidean 
lattice Z 4 . This system has very close nature with the gonihedric model because in both 



cases the action is proportional to the linear size of the manifold |L2) . In |T2| the authors 
where able to construct a transfer matrix for this theory as well. This transfer matrix 
describes the propagation of two-dimensional membrane M and our main result is that 
both theories can be solved exactly by using generalization of the Fourier transformation 
in loop and in membrane spaces. They represent two nontrivial theories in three and four 
dimensions respectively which can be solved exactly. 

In the next sections I shall present background material, the definition of the system 
of random surfaces in continuous space as well as on Euclidean lattice and describe its 
basic properties. In the third section the construction of the transfer matrix for the loops 
will be extended to the case of nonzero self-intersection coupling constant k. In the fourth 
main section we shall introduce the generalization of Fourier transformation in the loop 
space and demonstrate that loop Fourier transformation allows to diagonalize all transfer 
matrices which depend only on symmetric difference of loops. The net result is that all 
eigenvalues Ap of 3d transfer matrix are exactly equal to the loop correlation functions 
of the corresponding 2d system 

(^f^) 3d = < • • • ^ > 2d . (1) 

In this formula the loop momentum P is defined by the vectors f{, ...f* n (see Figure 8). In 
its form the last relation is very close with the ones discussed in ]Tj| [L7], |i~8|| . 

In the last section I shall consider the transfer matrix for membranes. This matrix 
also can be diagonalized by Fourier transformation and all its egenvalues are equal to the 
correlation functions of the corresponding 3d system. In particular we shall see that free 
energy of the membrane system is equal to the free energy of gonihedric system of loops 
and is equal to the free energy of 2d Ising model. 
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Figure 2: The random path and random surface with boundary points X and Y. In 
both cases the action is proportional to the linear size of the geometrical fluctuation 

2 Random Surfaces and Path Integral 

Surfaces in Continuous Space. Gonihedric string has been defined as a model of random 
surfaces with the action which is proportional to the linear size of the surface || |(| 

™ E X v ■ K - «0"I C (<(d-2)/d, (2) 

<ij> 

where \j is the length of the edge < ij > of the triangulated surface M and is the 
dihedral angle between two neighboring triangles of M sharing a common edge < ij > 
(see Figure 1,2) f\. The action was defined for the self-intersecting surfaces as well |J, 
it accounts self-intersections of different orders by ascribing weights to self-intersections. 
These weights are proportional to the length of the intersection multiplied by the angular 
factor which is equal to the sum of all angular factors corresponding to dihedral angles in 
the intersection (see Figure 3) 

Km ^2 ^ij ' ( k ~ a ijf + ••• + I 71 " — ct'ij^ )• (3) 

<ij> 

The number of terms inside parentheses depends on the order of the intersection r, where 
r is the number of triangles sharing the edge < ij >. Coupling constant k is called 
self- intersection coupling constant 0. Thus the total action is equal to 

A(M)=m A,, • - - <\,j - +mK ^ ■ ( |vr - aj/ + ... + |vr - ). (4) 

<ij> <ij> 

3 The angular factor defines the rigidity of the random surfaces and for £ < (d — 2)/d it increases 
sufficiently fast near angles a = tt to suppress transverse fluctuations [T(| . 
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Figure 3: The figure depicted the self-intersection of three planes and the associated 
action. In the given case a + (3 + 7 = 2ir. 

Both terms in the action have the same dimension and the same geometrical nature, the 
action (Q) measures two-dimensional surfaces in terms of length, self-intersections are also 
measured in terms of length (|3|). If k is small then self-intersections are permitted, if 
k increases then there is strong repulsion along the curves of self-intersections, and the 
surfaces tend to be self-avoiding 0. 

The model is a natural extension of the Feynman path integral for point-like relativis- 
ts particle (see Figure 1,2). This can be seen from (|2|) in the limit when the surface 
degenerates into a single world line ||Q 

Y Kj ■ k - «| ? < — > Y 

<ij> <ij> 

At the classical level the string tension is equal to zero and quarks viewed as open ends 
of the surface are propagating freely without interaction a c i assica i = 0. This is because 
the action (|2]) is equal to the perimeter of the flat Wilson loop A(M) —>■ m(R + T). 
As it was demonstrated in ||, quantum fluctuations generate nonzero string tension 
& quantum = (l — ^4), where d is the dimension of the spacetime, (3 is the coupling 
constant, a is a scaling parameter and <^ = in (||). In the scaling limit (3 — > (3 C = dj 'e 
the string tension has a finite limit. 

We shall also consider similar model of random membranes with the action which is 
proportional to the linear size of the corresponding worldvolume. This action for three- 
dimensional manifold has been defined as 

A{N) = m V \^-"ij\ C (5) 

<i,3> 

where A^- is the length of the edge < ij > of three-dimensional manifold iV which is 
constructed by gluing together three-dimensional tetrahedra through their triangular faces 
and 27r — uiij is the deficit angle on the edge < ij >, here Uij is the sum of the dihedral 
angles between triangular faces of tetrahedra which have the common edge < ij >. 

4 This property of the gonihedric action guarantees that spike instability does not appear here because 
the action is proportional to the total length of the spikes and suppresses the corresponding fluctuations 
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Equivalent spin system on a lattice j^j. In addition to the formulation of the theory 
in the continuum space the system allows an equivalent representation on Euclidean lat- 
tices where a surface is associated with a collection of plaquettes |7|]. The corresponding 
statistical weights are proportional to the total number of non-flat edges n^M) of the 
surface M as it is defined in (0) J7|. The edges are non-flat when the dihedral angle 
between plaquettes is equal to ir/2. The surfaces may have self-intersections in the form 
of four plaquettes intersecting on a link. The weight associated with self-intersections is 
proportional to 4mt-4 where n 4 (M) is the number of edges with four intersecting plaque- 
ttes, and k is the self-intersection coupling constant @ . The edges of the surface with 
self-intersections comprise the singular part of the surface. The edges of the surface where 
only two plaquettes are intersecting comprise the regular part of the surface. 

The partition function can be represented as a sum over two-dimensional surfaces M 
of the type described above, embedded into a three-dimensional lattice: 

Z(P)= e- 2 ? A{ - M \ (6) 

{surfaces M} 

where A(M) is the energy of the surface M constructed from plaquettes 

A gonihedHc (M) = n 2 (M) + 4«n 4 (M). (7) 

To study statistical and scaling properties of the system one can directly simulate 
surfaces by gluing together plaquettes following the rules described above, or to use the 
duality between random surfaces and spin systems on Euclidean lattice. This duality 
allows to study equivalent spin system with specially adjusted interaction between spins 
M. In three dimensions the corresponding Hamiltonian is equal to |7|] 

H gonihedric = 2^ ^f^f+S + ^ L a r a r+a+0 ^ 0r ^° r ^+a <T r+a+/3 <T r+ / 3' (°) 

where f is a three-dimensional vector on the lattice Z 3 , the components of which are 
integer and a, (3 are unit vectors parallel to the axes. The interface energy of this lattice 
spin system coincides with the linear action ([|), (^). 

The form of the Hamiltonian H(k) and the symmetry of the system essentially depend 
on k: when k ^ one can flip the spins on arbitrary parallel layers and thus the degeneracy 
of the ground state is equal to 3 ■ 2 , where iV 3 is the size of the lattice. When k = the 
system has even higher symmetry, all states, including the ground state are exponentially 
degenerate. This degeneracy of the ground state is equal to 2 3N [12| . This is because now 
one can flip the spins on arbitrary layers, even on intersecting ones. The corresponding 
Hamiltonian contains only exotic four-spin interaction term. 

The understanding of critical behaviour of this system comes from the analysis of the 
transfer matrix |]12| which describes the propagation of the closed loops Q in time direction 
with an amplitude which is proportional to the sum of the length l(Q) of the loop and of 
its total curvature k(Q). Because the amplitude essentially depends on length, as it takes 
place in 2d Ising model, this result directs to the idea that the system in 3d will show 
the phase transition which should be of the same nature as it is in 2d Ising system and 
that the transfer matrix describes the propagation of an almost free string of Ising model 



12]. Direct numerical Monte-Carlo simulations [19| and analytical solution of different 
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dim Q = dim M + dim E - dim R 



= 2 + (d-1)-d = 1 



Figure 4: The geometrical theorem provides an equivalent representation of the action 
A(M) in terms of total curvature k(E) of the curves Qi(E), Qk{E) which appear in the 
intersection of the (d-l)-dimensional plane E with the given two-dimensional surface M. 
The curvature k(E) should then be integrated over all planes E intersecting the surface 
M. 

approximations to the full transfer matrix confirm strong connection with the 2d Ising 
model 0. 

Below I shall remind the derivation of the corresponding transfer matrix for the case 



k — 0, which is based on geometrical theorem proven in |l I J and then I shall derive the 
transfer matrix for more general case k ^ 0. This will allow to study physical picture of 
string propagation which follows from the transfer matrix approach. 



3 Loop Transfer Matrix 

Geometrical Theorem In this section I shall review the construction of the transfer 



matrix for the system @ with k — | 12]| 



H gonihedric( K ~ °) ~ ~^ a r^r+a^ f + S+0 <7 ( 9 ) 

and then extend the result to nonzero case k ^ 0. 

In order to find transfer matrix for this system we have to use geometrical theorem 
proven in ||11|| . The geometrical theorem provides an equivalent representation of the 
action A(M) in terms of total curvature k(E) of the polygons which appear in the inter- 
section of the two-dimensional plane E with the given two-dimensional surface M (see 
Figure 4,5). This curvature k(E) should then be integrated over all planes E intersecting 
the surface M. 

The intersection of the two-dimensional plane E with the polyhedral surface M is the 
union of the polygons (see Figure 5) 

Q 1 {E),...,Q k (E) 
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k(E) = E | k - a E | 
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Figure 5: Intersection of the polyhedral surface M by the plane E. The image of the 
surface on the plane is a polygon P and its curvature is equal to k(E). 

and the absolute total curvature k(E) of these polygons is equal to 

k 

k(E)=Y,km= E k-«5l, (10) 

i=l <i,j> 

where afj are the angles of these polygons. They are defined as the angles in the inter- 
section of the two-dimensional plane E with the edge < ij > |Tl] (see Figure 5). The 
meaning of ( |T0D is that it measures the total revolution of the tangent vectors to polygons 
Qi(E), Qk(E). Integrating the total curvature k(E) (PH) over all intersecting planes E 
we shall get the action A(M) 



A(M) = — [ k(E)dE. (11) 
2tt J{e} 



Geometrical Theorem on a lattice ]7I /. One can find the same representation ([LI] ) 
for the action A(M) on a cubic lattice Z 3 introducing a set of planes {E x }, {E y }, {E z } 
perpendicular to x, y, z axis on the dual lattice Z\. These planes will intersect a given 
surface M and on each of these planes we shall have an image of the surface M. Every such 
image is represented as a collection of closed polygons Q(E) appearing in the intersection 
of the plane with surface M (see Figure 6). The energy of the surface M is equal now to 
the sum of the total curvature k(E) of all these polygons on different planes 

A(M)= £ k(E). (12) 

{E X ,Ey,E Z } 

The total curvature k(E) is simply the total number of polygon right angles. In the case 
when the self-intersection coupling constant k = is equal to zero we must compute the 
total curvature k(E) of these polygons ignoring angles at the self- intersection points! 
With (|1^) the partition function of the system can be written in the form 



Z((3) = Y,exp{-2f3j2HE)} (13) 

{M} {E} 
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where the sum in the exponent can be represented as a product 



exp{-2^k(E)} = IJe-W*) = J| e -W{E.) e -2/J*(^) J] e -2/W. (14) 

{£} {£} {iB,} {E y } {E x } 

The goal is to express the energy functional (0) and the product in terms of quantities 
defined only on two-dimensional planes in one fixed direction, let's say through {E z }. 

Transfer Matrix for k = j7I|/ . The question is: what kind of information do we 
need to know on these planes {E z } to recover the values of the total curvature k(E x ) and 
k(E y ) on the planes {E y } and {E x } for the given surface Ml 

Let us consider the sequence of two planes E\ and E l z +1 and denote by Qi the polygon- 
image of the surface M on the plane E\ and by Qi+i the polygon-image of M on the plane 
E\ +1 . The contribution to the curvature k(E x ) + k(E y ) of the polygons which are on the 
perpendicular planes between E\ and E\ +1 is equal to the length of the polygons Qi and 
Qi+i without length of the common bonds |L2] 



KQi) + KQi+i) - 2 ■ l(Qi n Q i+1 ) = l(Qi A Q i+1 ), (15) 

where the polygon-loop Qi A Q 2 = Qi U Q2 \ Qi H Q2 is a union of links Qi U Q2 
without common links Q\ D Q2- Therefore (see formula (p~2|) ) 

A(M) = ]T KE) = E + l ^ A Qi+i). (16) 

{E} {E z } 

The partition function (|H|) can be now represented in the form 

Z{p)= E Kp(Q 1 ,Q 2 ).-.K fi (Q N ,Q 1 )=trKjf, (17) 
{Qi,Q 2 ,-,Qiv} 

where Kp{Q\, Q 2 ) is the transfer matrix of size 7x7, defined as 

K 9 K To hednc (Qi, Q2) = exp{-(3 [k(Qi) + 2/(Q x A Q 2 ) + k(Q 2 )] }, (18) 

where Q\ and Q 2 are closed polygon-loops on a two-dimensional lattice T 2 of size N x N 
and 7 = 2 N2 is the total number of polygon-loops on a lattice T 2 . The transfer matrix 
( |I8"D can be viewed as describing the propagation of the polygon-loop Qi at time r to 
another polygon-loop Q 2 at the time r + 1 (see Figure 7)|. 

The functional fc(Q) is the total curvature of the polygon-loop Q which is equal to the 
number of corners of the polygon ( the vertices with self-intersection are not counted!) 
and l(Q) is the length of Q which is equal to the number of its links. The length functional 
l(Qi AQ2) was defined above (|T5|) . The operation A maps two polygon-loops Qi and Q 2 
into a polygon-loop Q = Qi A Q 2 0. The length functional 1{Q\ A Q 2 ) defines a distance 
between two polygon-loops Qi and Q 2 . 

5 Layer-to-layer transfer matrices for three-dimensional statistical systems have been considered in the 
literature |2^|. Using Yang-Baxter and Tetrahedron equations one can compute the spectrum of 
the transfer matrix in a number of cases pl| . In the given case the transfer matrix has geometrical 
interpretation which helps to compute the spectrum. 

6 Note that the operations U and n do not have this property. These operations acting on a polygon- 
loops can produce link configurations which do not belong to II. The symmetric difference of sets Qi AQ2 
is an important concept in functional analysis p3j. 
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k(E)=16 k(E)=20 k(E)=20 

z y x 

Figure 6: The set of planes {E x }, {E y }, {E z } perpendicular to x,y,z axis intersect a 
given surface M in the middle of the links. On each of these planes we shall have an 
image of the surface M. Every such image is represented by closed polygon-loops P(E) 
appearing in the intersection of the plane with surface M. The energy of the surface M 
is equal to the sum of the total curvature k(E) on all these planes as it is given by the 
formula (|L2D . For the surface on the picture this sum is equal to A = 16 + 20 + 20 = 56 
times a|. 



P 2 



I PAP 2 

Figure 7: The transfer martix (fL8l) can be viewed as describing the propagation of the 
polygon-loop P\ at time r to another polygon-loop P2 at the time r+ 1. For the surface on 
the picture we have k{P\) = A,k(P 2 ) = 12,/(0APi) = 12,/(PiAP 2 ) = 16,/(P 2 A0) = 12. 
Summing all these quantities in accordance with the formula Qlfip we shall get A(M) = 
4 + 12 + 12 + 16 + 12 = 56 times a|. This, as it should, coincides with the previous result. 
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The eigenvalues of the transfer matrix K/3(Qi,Q 2 ) define all statistical properties of 
the system and can be found as a solution of the following integral equation in the loop 
space II [[TTJ 

E k p(Qi> Q*) ^(Q2) = A(/3) *(Q0, (19) 
{q 2 } 

where ^(Q) is a function on loop space. The Hilbert space of complex functions ^(Q) on 
n will be denoted as H = L 2 (U). The eigenvalues define the partition function ( |T3"D 

Z* d ((3)=A» + ...+A%, (20) 

and in the thermodynamical limit the free energy is equal to 

" P MP) = & ^ ^ Z 3d (P) = Jim ± In A . (21) 

Finite time propagation amplitude of an initial loop Qi to a final loop Qf for the time 
interval t = Mj [3 can be defined as 

K(Q l ,Q f )=A M Yl Kp{Q h Q x ).--Kp{Q M ^Q f l (22) 

{Qi,Q2,---,Qm-i} 

where we have introduced natural normalization to the biggest eigenvalue A and M < N . 

Loop Transfer Matrix for k / 0. We would like to generalize the construction of 
the transfer matrix to the case of nonzero self-intersection coupling constant k. As we 
already mentioned above surfaces on a lattice may have self-intersections and we associate 
the coupling constant k with these self-intersection edges of the surface. When we consider 
the sections of the lattice surface M by the planes E the self-intersection edges will appear 
on a plane as polygon-loops with self-intersections. When k = we have been instructed 
to compute the total curvature k(E) of these polygon-loops simply ignoring angles at the 
self-intersection points. Now when k ^ we have to take into account self-intersection 
vertices of polygon-loops. 

The contribution to the action which comes from simple polygon right angles is the 
same as before ( [TED 

J2 k{Qi) + l{QiAQM). (23) 

{Em} 

The contribution coming from self-intersections should be divided into contribution com- 
ing from vertical and horizontal edges of the surface. The vertical and horizontal edges of 
the surface are defined relative to the "time" planes E z . The contribution from vertical 
edges is simply the number of self-intersection vertices of loops on E z planes 

{Ez} 

where k int {Q) is the number of self-intersection vertices of loop Q. 

To count contribution from horizontal edges of the surface one should introduce the 
orientation of the polygon-loops on every E z plane. The same orientation is ascribed 
to all polygon-loops on a given plane E z : clock-wise or anti-clock- wise. Orientation is 

automatic if Q is thought as interface between spins. Conjugation Q — ► Q is equivalent 

to Z 2 spin flip. So there is Q and Q with opposite orientation. 
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Now let us consider oriented polygon-loops on two planes E\ and E l z +1 . Then one 
should select that bonds on these polygon-loops which are common. Only bonds with 
the opposite orientation should be counted. The number of common bonds with opposite 
orientation we shall denote as l{Qi fl Qj). The notation H is used to denote the common 
bonds, the symbol <-> on the top of it to denote that they should have opposite orientations. 
Therefore the contribution from horizontal edges can be expressed as 

Ak ■ ]T l(Qi n Qj) (25) 

{E z } 

and the contribution from all self-intersections as 

4«:- £ k int (Qi) + KQiriQ i+1 ). (26) 

The total action (0) now takes the form: 

A(M) = KQi) + KQi A Qi+i) + 4k • [ ^(QO + n Q i+1 ) ]. (27) 

The partition function (|IJ|) can be now represented in the same form (|I7|) where K@(Qi, Q 2 ) 
is again the transfer matrix of size 7x7, defined as 

Kp(Q 1 ,Q 2 ) = exp{-(3 [k(Qi) + 2l(Q 1 A Q 2 ) + k(Q 2 )]} ■ 

exp{-4K(3 [k int {Qi) + 2/(Qi n Q 2 ) + fc int (Q 2 )]}, (28) 

where Q\ and Q2 are closed oriented polygon loops on a two-dimensional lattice T 2 of size 
N x N and 7 is the total number of polygon- loops on a lattice T 2 . Let us now consider 
the limit 

K — ► OO, 

then from our expression for the transfer matrix (|2q) one can conclude that on every time 
slice there will be no polygon-loops with self-intersections. This means that they are just 
nonsingular self-avoiding oriented loops. In addition there is strong correlation between 
loops on consecutive planes: common bonds with opposite orientation are also strongly 
suppressed. Which means that propagation of oriented loop in time has small probability 
to invert its orientation. 



4 Loop Fourier transformation 



In the approximation when we drop the curvature term in the expression (IB) for the loop 
transfer matrix K(Qi, Q 2 ) the spectrum of the matrix has been evaluated exactly in the 
articles 
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|12j. In this section I shall present an alternative method of derivation which 
is based on generalization of Fourier transformation and allows to carry out analogy with 
quantum mechanics of point particles and to emphasize generality of the approach. 

I shall consider two different approximations to the initial transfer matrix (1TB]) sug- 



gested in JT2| and [0. In the first one it was suggested to substruct the term 2k(Q\ UQ 2 ) 
and to consider the transfer matrix of the form flT2"| 



K{Q 1 A Q 2 ) = exp{-p [k(Qi A Q 2 ) + 2Z(Qi A Q 2 )\ }, 



(29) 
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and in the second one to consider the transfer matrix |p~3 | 



K{Q X A Q 2 ) = exp{-2(3 l(Q 1 A Q 2 ) }■ (30) 

In both cases the transfer matrix depends only on the symmetric difference Q\ A Q 2 and 
the integral equation (|T9|) take the form 

£ JWi A Q 2 ) ¥(Q 2 ) = A({3) iif(Q 1 ). (31) 
{Q2} 

The fact that transfer matrix now depends only on the symmetric difference is of paramount 
importance in solving these models. The reason is that we can make a change of variable, 



from Q 2 to Q at fixed Q\ [13 



Q = Q l A Q 2 (32) 
and see that equation transforms to the form 

J2 K P (Q) M>(Q A = A(/3) (33) 
{Q} 

The integral equation will be solved if we can find the set of loop functions *&{Q) which 
satisfy the following functional equation on loop space II 

*(QAQ 1 ) = *(Q)-*(Q 1 ). (34) 

For these functions the equation fl3"3|) reduces to 



y £K p {Q)*(Q)=A{p) (35) 
{Q} 

and we can find all eigenvalues. The constant function ^ = 1 is the solution of the 
functional equation (|3~4D and due to Frobenius- Perron theorem corresponds to the largest 
and non-degenerate eigenvalue A0 of the symmetric transfer matrix K(Qi A Q 2 ). 

The functional equation (0) is analogous to the equation ip(x+y) = i()(x)-if)(y) appear- 
ing in ordinary quantum mechanics, which allows to find all eigenvectors and eigenvalues 
of the integral equation 



Kp(x - y)^/j(y)dx = \(P)i[>(y), K p (x - y) = exp(-(3(x - y) ) (36) 
with the well known solutions 

rj)(x) = e ipx , X p (x) = e- p2/4f3 . 



Our idea is to find all solutions of the functional equation ( p4|) and the analog of the plane 
wave solution in loop space. From functional equation ( |34"D it follows that if Q = Qi then 



*(0) = ^ 2 (Q) (37) 
and if Q = then \I/(0) is equal to zero or to one, therefore nontrivial solution is 

* 2 (Q) = 1 ^(Q) = ±l. (38) 



14 



This suggests the following form of the eigenfunctions 

tfp(Q) = e i«(PnO) (39) 

which are numbered by the loop momentum P and take values ±1. The s(Q) is the area 
of the region with the boundary loop Q. These solutions are in pure analogy with plane 
waves in quantum mechanics . To check that they are indeed solutions of the functional 
equation ([34]) we have to use the set identity 

(Qi A Q 2 ) n P = (Qi n P) A (Q 2 n P) (40) 

which is easy to verify and the relation s(Qi) + s(Q i+ i) — 2 • s(Qi fl Qj+i) = s(Qi A Qi+i)- 
To carry out the analogy with plane wave solution of quantum mechanics deeper one 
should interpret the loop P as a loop momentum conjugate to a loop coordinate variable 
Q. To prove orthogonality of these loop wave functions let us compute the product 

. ^ Pa = £ * Pi (Q) ^ P2 (Q) = J2 e i ^ nQ)+i ^ nQ) = £ e ^« p i Ap 2)nQ) ( 41 ) 
«} Wl {Q} 

where we again used the above set identity and the relation (|TJ) for s(Q). The last sum 
is equal to zero if Pi A P 2 7^ and is equal to the number of loops 7 on the lattice if 
P 1 AP 2 = 0, thus 

e i-((PiAP 2 )nQ) = 2 n2 5(P 1 , P 2 ). (42) 

{<?} 

The last statement follows from the fact that the product between the constant function 
^0 = 1, corresponding to the largest and non-degenerate eigenvalue A0, and all others 
exp(iirs(P HQ)) is equal to zero if P is nonempty 

J2 1 . e i*-((PnQ) = 0, if P^ds. (43) 
{Q} 

The number of functions in this orthogonal set is equal therefore to the number of closed 
loops P, that is 7 = 2 N . 

Having in hand all eigenfuctions we can find all eigenvalues from ([35]) for transfer 
matrices © and © 

A p = e - i7TS ( pnQ) -P [k ( Q)+2im A P = J2 e - i7TS ( pnQ) - 2f3l( > Q \ (44) 

{Q} {Q} 

Because the expression Z){q} e~ 2)31 ^ is the partition function of the 2d Ising model 0, [| 
and J2{q} e~^ 2l ^ +k ^ is the partition function Z 2d of the two-dimensional model solved 
in [12j we see that the largest eigenvalue A0 is exactly equal to a corresponding partition 



function 

A = Z 2d A = Z 2d ~ Isin9 . (45) 

For the ratio of eigenvalues we shall have [] 

^£ _ e -i™(Pr\Q)-2pl(Q) i^d-Ising = < e -ins(PnQ) >q ^q) 

Ae {Q} 



7 In what follows we shall consider the transfer matrix (|3y) . Analogous formulas are valid for the 
transfer matrix (p9|) with curvature term k(Q). 
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which is identically equal to the spin correlation functions of the 2d Ising model 



A iV,...,n 

^ = < (Tfi ■ ■ ■ CT r - >2dlsing ■ (47) 

A 

In this formula the loop momentum P is defined by the vectors f{, ...?v In particular for 
the smallest (one-box) loop momentum P = □ we shall have 

= < a? > = h(J3) (48) 

which coincides with the magnetization of 2d Ising model. Because it does not depend 
on position of the loop on the lattice, this eigenvalue is iV 2 degenerate. For the two-box 
loop momentum P = consisting of two neighboring boxes we shall have 



Aw 



< a&rZx > = - U (P) ( 49 ) 



which is 2A^ 2 degenerate and after summation over all 2iV 2 bonds of the lattice we shall 
get the internal energy of the 2d Ising model. 

We can also use loop wave functions to define generalized loop-Fourier transformation 

as 

*(^) = iE eWS(PnQ) *W) (5°) 

1 {Q} 

and then compute the loop Fourier transformation of the transfer matrices (|29|) and (|3(i|) 

as: 

K(P 1 ,P 2 )= ^ K{Qi A Q 2 ) e i7rs ( P l n< 3l)-^(-P2nQ2) _ yy^ ^ e iirs(P 1 nQ2)-iirs(P 2 riQ2) 

{Qi,Qa} {Q2} 

= A Pl 8(P 1 ,P 2 ). (51) 

The finite time propagation amplitude of an initial loop Pi to the final loop Pf for the 
time interval t = Mj (3 ( p2"D we have 

K(P h P f ) = (^) M 5(P h P f ) = e Mln{ ^ ] 5(P h P f ) = e~ MH ^ 6(P h P f ), (52) 

therefore for the loop Hamiltonian we shall get 

H(P) = -ln(^) (53) 

where H(P) is the energy of the coherent state with loop momentum P. 

Spin system which corresponds to the matrix K(QiAQ i+ i) (|5^. The natural question 
which arises here is how nontrivial is the 3d system which has been solved in this way? To 
make the answer as clear as possible we should go back and express the system which has 
been solved again in terms of spins on a 3d lattice and to see what is left out of original 
interactions. 

We ignore the curvature terms k(Qi) in the transfer matrix, this means that we ignore 
the contribution coming from the edges of the surface which are along the time direction, 
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r r r r r 

1 2 3 4 5 



Figure 8: The loop momentum Pf\ r ..f n is defined by the vectors f{,...,r\ l and the 
corresponding eigenvalue Ap^ ^ is equal to the correlation function of the n spins 

< Gf\ ■ ■ ■ Gf* n >2d- 



or z direction. The contribution coming from the edges which are parallel to the x and y 



direction are left and their contribution is equal to (15) 



]T l(Qi AQ i+1 ). (54) 

{E z } 

In terms of spin variables the original action can be expressed as a sum of energies coming 
from all edges of the surface (see formulas (5a) and (7) in 0) 

H= J2 H ed 9 e, (55) 

over all edges 

where 

Hedge = U\U -\ + U X U-l = Ox0 2 -lO "_ 2 . (56) 

is a four-spin interaction term with spins distributed around the edge. Because we ignore 
the edges in the z direction, we have to sum in fl5"5|) only over A e d g es in x and y direction. 
This corresponds to the summation over four-spin interaction terms (|56|) which lie only 
on the planes E x and E y . These planes are normal to the x and y direction. Therefore our 
approximation corresponds to the 3d lattice spin system in which four-spin interactions 
take place only on vertical planes E x and E y 

H Qi aq 2 = (57) 

It is obvious that this spin system can not be factored into non interacting two-dimensional 
subsystems. We can ask what is the effective interaction between spins which lie on a 
given two-dimensional E z plane? Let us label spins which are on that plane as er, spins 
which are on the previous plane as A, and spins which are on the next plane as /i. The 
part of the Hamiltonian which contains only a spins can be written in the form: 

X] ^iOi+l^A+l + = + HiHi+lj&i&i+l = Jeff&i<Ti+l, (58) 
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where the effective coupling J e ff depends on spins on the neighboring planes. It takes the 
values —2,0,2 with probabilities 1/6,4/6, 1/6 simply because 

Jeff = W+l + f'if'i+l- 

For the observer who lives on the plane E z spin interactions look like 2d Ising model with 
randomly distributed coupling constant J e ff- 

3d Ising Transfer Matrix. Similar expression for the transfer matrix has been derived 
in [12j for the 3d Ising model 

Kf- Isin \P u P 2 ) = exp{-(3 [Z(Pi) + 2s(P 1 A P 2 ) + l(P 2 )) }, (59) 

where s(P) is the area functional. It is instructive to compare 3d-gonihedric and 3d Ising 
systems. One can see simple hierarchical structure of the corresponding transfer matrices 
( |I~8D , (0). If one orders the functionals in the increasing geometrical complexity 

curvature k(P), length l(P), area s(P), (60) 

then it becomes clear that gonihedric transfer matrix depends on the first two functionals: 
curvature and length, while the 3d Ising model depends on the last two functionals: length 
and area. Thus if one puts the models in order of increasing complexity one should 



consider gonihedric model and then Ising model [12 



Let us also consider similar to (|29|), (|30|) approximations to the 3d Ising transfer 
matrix, they are: 

K{Q 1 AQ 2 ) = ex P {-f3 [l(Q 1 AQ 2 )+2s(Q 1 AQ 2 )] }, K(QiAQ 2 ) = exp{-(3 2s{Q l AQ 2 ) }, 

(61) 

Both approximations can be solved by using the set of eigenfunctions (|39| ) 



A p = V" e -iTrs(PnQ)-i3[i(Q)+MQ)] a p = e~ ins ( PnQ) ~ 2f3s{Q) . (62) 

{Q} {Q} 

The first relation expresses the eigenvalues through the 2d Ising model in background 
magnetic field, the model which has not been solved yet exactly |HJ, and in the second 
case through free spins in background magnetic field. 



5 Matrices depending only on symmetric difference 

It is obvious that the main reason why we have been able to reduce solution of three- 
dimensional statistical system to a solution of two-dimensional system is special depen- 
dence of the transfer matrix argument from symmetric difference Q\ AQ 2 . 

We are aimed to formulate the following general result: all transfer matrices of three- 
dimensional statistical systems which depend only on symmetric difference Qi A Q 2 

Z 3d (f3) = TrK N = e~^^ N \ K{Q 1 A Q 2 ) = e -Wi**0 (63) 

where Q is arbitrary energy functional, can be diagonalized by using orthogonal set of 
functions ^p(Q) = e l7Ts( - Pn ®\ The eigenvalues can be found exactly in the same way as 
we described above and are equal to the following averages 

A P = e~ i7Ts{PnQ) K(Q) = J2e~ lns{PnQ) ~^ Q) . (64) 
{Q} {<?} 
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The eigenvalues Ap of three-dimensional system are exactly equal to the correlation func- 
tions of the two-dimensional system with the partition function 

Z 2d ({3) = ]T = e-PfM-^ = \q + ... + A 2 V X , (65) 

{Q} 

where fzdift) is the free energy of the 2d model and A, are the eigenvalues of the corre- 
sponding transfer matrix. In particular from flB"3]) we see that the largest eigenvalue A0 is 
equal to the partition function of 2d system 

A = E^ (Q) = Z 2d (f3) = A^ + ... + A^_ 1; (66) 
{<?} 

therefore we can express free energy of 3d system in terms of largest eigenvalue A0 and 
then through the partition function of the corresponding two-dimensional system 

= Jim -L ln{[\$ + ... + A^J" + ...} = -13 MP). (67) 
Thus we have amazing equality 

which proves identical critical behaviour of 2d and 3d systems! Because 

^ e - l ns(Pn Q )-f3n(Q) /Z 2 d{(3) = < e -i„s(PnQ) >q= = <(Tfi ... aK >M; 
{Q} 

from (|4]) it follows that the ration of eigenvalues is equal to: 

= <**..• <r P - > 2rf . (68) 

A0 

In this formula the loop momentum P is defined by the vectors f[, ...7%,. It is clear that the 
correlation functions of the three-dimensional system are more complicated conglomerates. 

It is easy to understand that the result is even more general: any d — dimensional 
system which has transfer matrix depending only on symmetric difference Q1AQ2 can be 
reduced to a corresponding (d — 1) — dimensional system and its eigenvalues are simply 
the correlation functions of the lower dimensional system. Below we shall consider an 
important example of such dimensional reduction in four-dimensions. 



Transfer Matrix for Membranes In [12| the authors have constructed the transfer 



matrix in four-dimensions which describes a propagation of two-dimensional membrane. 
Below we shall see that this system also can be solved in terms of three-dimensional 
system considered in previous sections. 

Transfer matrix which describes the propagation of two-dimensional membrane M has 
the following form |L2| 



KpiMx, M 2 ) = exp{-(3 fr(Mi) + 2A{M 1 A M 2 ) + X (M 2 )] }, (69) 
where x(M) defines the module of the Euler character of the surface M 

X (M) = ]T|27r-^|, (70) 
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summation is over all vertices of triangulated surface M and A(M) is the gonihedric 
functional which we defined earlier. If we introduce the set of orthogonal functions 

* P (M) = e lnv{PnM \ (71) 

where v(M) is the volume of the region whose boundary is the surface M, then transfer 
matrices which depend only on symmetric difference of surfaces M\ and M 2 

K(M 1 ,M 2 ) = exp{-p [ X (M 1 AM 2 )+2A(M 1 AM 2 )}}, K(M 1 ,M 2 ) = exp{— 2/3 A(M 1 AM 2 )}, 

(72) 

can be diagonalized by using the above set of functions ([H|) 

A p _ e -iTvv(PnM)-f3[ x (M)+2A(M)] ^ \ p = ^ e -wv(PnM)-2/3A(M)_ (73) 
{M} {M} 

For P = we have 

A = £ e -W(M) (74) 

{M} 

and it coincides with the partition function of the gonihedric system (|J). The free energy 
of this four-dimensional model is equal therefore to 

UM = fsM = f2d-Isin 9 {P) (75) 

and the eigenvalues (|73|) are the correlation function of the 3d gonihedric system 

= {Correlation functions) . (76) 



A y 

Thus the eigenvalues of the transfer matrix in four-dimensions can be described by the 
correlation functions of the system of one less dimension. 



6 Discussion 

We have seen that transfer matrix approach allows to solve large class of statistical systems 
when transfer matrix depends on symmetric difference of propagating loops, membranes 
or p-branes. The eigenvalues are equal to the correlation functions of the corresponding 
statistical system, which has smaller dimension than the original one. This dimensional 
reduction expresses the hierarchical structure in which more complicated high-dimensional 
systems are a superposition of lower-dimensional ones WQ. This structure is especially 
visible when one uses generalization of Fourier transformation in loop space or in p-brane 
spaces. 

In conclusion I would like to acknowledge R.Kirschner, B.Geyer and W.Janke for 
discussions and kind hospitality at Leipzig University. This work was supported in part 
by the EEC Grant no. HPMF-CT-1999-00162. 
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